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We calculate correction to the nuclear fission barrier produced by the atomic electrons. The result 
presented in analytical form is convenient to use in future nuclear calculations. The atomic electrons 
have a small stabilizing effect on nuclei, increasing lifetime in nuclear fission channel. This effect 
gives a new instrument to study the fission process. 
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PACS numbers: 24.75.+i,31.30.Gs 

In this paper we consider the effect of atomic electrons 
on nuclear fission. This effect is small for all known nu- 
clei. However, it might still be detectable. It can be ob- 
served by comparing fission probabilities of nuclei striped 
from all electrons with those in neutral atoms. This may 
reveal important information about fission. The effect is 
expected to play more significant role in fission of super- 
heavy nuclei leading to increased lifetime of the nuclei. 
The study of the superheavy elements is now a popular 
area of research due to the search for the hypothetical 
stability island (Z=114 to 126) and good progress in the 
synthesis of the elements on accelerators up to Z = 118 
(see, e.g. reviews It is natural to expect that 

physics of the superheavy nuclei has some phenomena 
which never manifest themselves for lighter nuclei. 

In this paper we perform relativistic atomic many- 
body calculations and give a formula which fits calculated 
total energy of an atom as a function of nuclear charge 
Z and nuclear radius Rn- This formula can be used in 
nuclear calculations of the fission probabilities. 

Nuclear fission is the main decay channel for heavy nu- 
clei starting from H^'Cf iQ]. During the fission a nucleus 
goes throw several stages of deformation before separa- 
tion of the fragments. The process is accompanied by the 
nuclear rotation. 

The electrons are far away from the nucleus and the 
electron density does not rotate when the nucleus rotate. 
Thus, the electrons feel nuclear charge density averaged 
over the nuclear rotation. This nuclear charge density 
for spinless nucleus is spherical. For a non-zero nuclear 
spin there is an electric quadrupole moment correction 
which is negligible anyway. Here there is the difference 
with the muon case considered in Ref. 0]: the muon 
is inside the nucleus, so the calculation in the rotating 
(frozen nucleus) frame is appropriate. 

Therefore, the change of the electron energy during the 
fission process is similar to the well known volume (or 
field) isotope shift of atomic transition frequencies which 
is determined by a change of a single nuclear parameter, 
mean squared nuclear charge radius. Total electron en- 
ergy is negative and it goes up when nuclear radius is 
increasing. This means that electrons take some energy 
from the process making fission more difficult. 

Stress that in this paper we do not perform any nu- 
clear calculations. A real shape of the fission barrier is 
complicated, with several minima, and an accurate cal- 




FIG. 1: Nuclear energy U{a) as a function of deformation 
parameter a. 



culation would require sophisticated nuclear codes. How- 
ever, to explain the physics involved and find out if the 
measurements are feasible we start from presenting some 
rough estimates using a simplest parabolic barrier model 
(Fig. [J) [1, 0| ■ The probability of the tunneling through 
the barrier is given by 0] 



P = 



1 + exp I 



huj > 



(1) 



where C/s is the maximum of the potential energy, (point 
B at Fig.p and is of the order from 0.5 to I MeV 0. 
For the spontaneous fission the tunneling amplitude is 
very small. Indeed, according to ref. [6| C/b — Va ~ 
5MeV. Therefore, 2-k\U a - VsM'^uJ > f and 



P « exp 27r 



E -Ub 



(2) 



for energies not too close to Ub ■ To estimate the effect of 
electrons on nuclear fission we should take into account 
the electron energy Ee{a) = Ee{aA) + k{a — a^) + ... 
where a is the deformation parameter (see Fig. [1]). The 
constant Ee{aA) does not change the probability of the 
tunneling since the constant potential does not influence 
the wave function (the bound state energy E and the 
potential U change by the same amount Ee{aA) so the 
difference E — Ub does not change). The linear term 
k{a — aA) can be incorporated into the effective parabolic 
potential and slightly changes the difference E — Ub- In- 
deed, it vanishes at the point of minimum a ~ aA where 
the bound state is located, and modify the potential near 
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the maximum which determines the tunnehng probabil- 
ity. A simple calculation gives the following result for the 
tunneling probability 



5E 

P = Poexp( -27r— 



(3) 



where 5E = Ee{ctB) — Ee{aA) is the change of the elec- 
tron energy. This change is positive and therefore the 
factor in ^ is smaller than one. This means that the 
electrons make the probability of nuclear fission smaller. 

To calculate 5E we need to calculate total electron en- 
ergy of an atom for two different nuclear charge distribu- 
tions corresponding to point A and point B (Fig. [T]). In 
light actinide atoms, e.g. ^^^J7, the deformation changes 
from 0.2 in the minimum to about 1 in the last maximum 
of the fission barrier (see e.g. @,@|). In heavier elements 
the difference is significantly smaller, about 0.2-0.4. As 
we mentioned above, after the averaging over nuclear ro- 
tation the charge distribution is spherically symmetric. 
For an estimate we assume that the change of the radius 
of the nuclear charge distribution from point A to point B 
is 10% (an accurate calculation will be discussed below). 
We calculate nuclear potential by integrating standard 
Fermi distribution for nuclear density 



p{r) 



C 



1 



exp 



. D , 



(4) 



where D = d/41n3, d=2.3 fm, ~ 1.2(3Z)i/3 fm and 
C is a normalization factor defined by / pdV = Z\e\. 

To find the electron energy of the atom we solve self- 
consistently a set of relativistic Hartree-Fock equations 
for single-electron orbitals (atomic units) 



+ -Mr)~ 2 + a2(e,;-F) <?.(r) = 0, 



dr r 

-7- - -Mr) + (e 



dr 



V)h{r) = 0, 



(5) 



here n = (— l)'^-'^^^^(i + 1/2),/ and j are angular and 
total electron momenta and V is the sum of nuclear po- 
tential, found by integration of nuclear density Q and 
the self-consistent Hartree-Fock potential of atomic elec- 
trons. Index i numerates single-electron states. 

Total electron energy in first-order in Coulomb inter- 
action is given by 



z 



=1 



(6) 



where are Hartree-Fock energies (Eq. ([5|)) of Z atomic 
electrons, q{ijij) = — qiijji) and q[abcd) is a 

Coulomb integral 

q{abcd) =11 i'a{ri)ipb{r2)-. — rV'c(f'i)V'd(»^2)c?i"irfr2- 

J J Fi - r2 

The first-order Coulomb correction is important because 
it excludes double counting of the energy of the Coulomb 



interaction between electrons. Indeed, the energy of the 
Coulomb interaction between electrons i and j is included 
in both €i and e^. 

The results for Z from 80 to 160 are presented in Ta- 
ble|Tl In this table we present the change of the total elec- 
tron energy of the atom and corresponding change in the 
fission probability when nuclear radius changes by 10%. 
Corrections to the fission probabilities were calculated at 
TiLj = 0.5 MeV using formula Q. At = 1 MeV the 
change of the probability is about 2 times smaller as is 
obvious from Eq. Eqs. ^ and ^ start to devi- 

ate significantly from each other at about Z = 130. In 
the vicinity of the stability island (Z « 120 to 130) elec- 
trons decrease fission probability by few per cent. The 
effect is large at extremely high Z leading to doubling of 
the nuclear lifetime a,t Z = 160. The Is electrons give 
about 30% of the total energy and more than a half of 
its change due to the change of the nuclear radius. This 
is important because the total energy depends on the 
electron configuration and accurate calculations should 
include investigation which configuration corresponds to 
the ground state. Large contribution from the Is elec- 
trons show that the effect of external electrons can be 
neglected. 

Since the change of the electron energy is dominated 
by Is electrons it is instructive to consider a simple pic- 
ture with only one electron in the Is state. The change of 
its energy due to the change of the nuclear radius is very 
small. Therefore it is natural to try perturbation the- 
ory. First order correction to the Coulomb energy can be 
written as 



AE{rN) = -eJi 



- Ze/r)%l?{r)dV, 



(7) 



where (j) is the nuclear potential corresponding to the 
finite nuclear radius r^v. The change of the energy due 
to small change of nuclear radius can be found as the 
difference 



5{AE) = AE{rNB) - AE(rNA) 



(8) 



Here ttva and r^^B are nuclear radii at points A and B 
of Fig. [TJ The results for the relative change of the ener- 
gies of Is electrons for different Z obtained with the use 
of the relativistic Coulomb wave functions in the Eq. ([7]) 
are presented in second column of Table [Hi The change 
of nuclear radius is 10% and the radii are the same as in 
Table [D The third column of Table |TT] presents similar 
results which were obtained by replacing the Coulomb 
wave functions with the numerical solutions of the Dirac 
equation corresponding to finite nuclear radius. The re- 
sults are very close for small Z but deviate significantly 
for higher Z . This actually means that the perturbation 
theory doesn't work in spite of the fact that the correction 
to the energy is very small. The reason for this is very 
simple: the perturbation is not small but it is localized 
in the very small volume of space limited by the nuclear 
radius. Indeed, only the nuclear volume contributes to 
the integral ([7]). The finite nuclear radius dramatically 



3 



TABLE I: Energy of Is electron and total electron energy (MeV) of many-electron atoms as a function of nuclear charge Z and 
nuclear radius rjv. In the last column we present rough estimates of the relative change of fission probability assuming 10% 
change of nuclear radius between points A and B; the actual increase with nuclear charge is slower since the relative change of 
radius in heavy elements is smaller than in light elements. 



z 


TjSf 


Energy [MeV] 


TjSf 


Energy [MeV] 




^tot 






[fm] 


Is 


Total 


[fm] 


Is 


Total 


[MeV] 


[MeV] 


[%] 


80 


6.44 


-0.083656 


-0.534682 


7.09 


-0.083649 


-0.534664 


0.00001 


0.00002 


0.02 


90 


7.75 


-0.110411 


-0.721344 


8.53 


-0.110387 


-0.721283 


0.00002 


0.00006 


0.08 


100 


8.03 


-0.142969 


-0.948529 


8.84 


-0.142908 


-0.948367 


0.00006 


0.00016 


0.20 


110 


8.29 


-0.183076 


-1.225050 


9.12 


-0.182917 


-1.224613 


0.00016 


0.00044 


0.55 


120 


8.54 


-0.233429 


-1.563897 


9.39 


-0.233024 


-1.562712 


0.00041 


0.00118 


1.5 


130 


8.76 


-0.298439 


-1.987579 


9.64 


-0.297409 


-1.984307 


0.00103 


0.00327 


4.0 


140 


8.99 


-0.385348 


-2.538006 


9.89 


-0.382769 


-2.528802 


0.00258 


0.00920 


11 


150 


9.20 


-0.505625 


-3.289468 


10.12 


-0.499261 


-3.263392 


0.00636 


0.02608 


28 


160 


9.40 


-0.675140 


-4.345553 


10.34 


-0.660488 


-4.278862 


0.01465 


0.06669 


57 



TABLE IL Relative change (in per cents) of the energy of the 
Is state due to change of nuclear radius by 10%; (1) integra- 
tion of Eq. ([7} with Coulomb Is functions; (2) integration of 
Eq. ((TJ with the functions obtained by numerical solution of 
the Dirac equation with finite size nucleus; (3) the difference 
between solutions of the Dirac equation for two nuclear radii; 
(4) formula from Eq. 



z 


(1) 


(2) 


(3) 


(4) 


80 


0.01 


0.01 


0.01 


0.01 


90 


0.03 


0.02 


0.02 


0.02 


100 


0.07 


0.04 


0.04 


0.04 


110 


0.16 


0.08 


0.08 


0.09 


120 


0.47 


0.17 


0.16 


0.23 


130 


1.74 


0.35 


0.32 


0.82 



changes the electron wave function inside this volume. 
In the end we use accurate numerical calculations rather 
than perturbation theory to find the change of electron 
energy due to change of nuclear radius. 

In fourth column of Table [ll] we present the results 
of the numerical solution of the Dirac equation with the 
finite nuclear size. These results are in good agreement 
with the data in previous column. The energy shift is 
dominated by the s-wave electrons. Note that the single- 
electron s-wave energy shift due to finite nuclear size can 
be approximated by a semi-empirical formula which is 
accurate to few per cent from Z = 1 up to Z = 100: 



V 57(7 + 1) 



2Z 



ao 



27 



(9) 



where v is the effective principal quantum number {v = 
Zeff^a.u./ {2\E^\); Z^ff — Z, v — 1 ior the Is state in 
the Coulomb field; Z^ff — I, ly ^ 1.7 for an external 
electron in a neutral atom), 7 — y^l — (Za)^, uq is Bohr 
radius. If change of radius is small {Sr^ ^ fjv) then 
Eq. © leads to the following expression for the volume 



isotope shift 
6{AE,) 



\E,. 



V 5(7+1) 



2Z 



(10) 



The energy shifts calculated with the use of Eqs. ([9]) and 
PU)) are presented in last column of Table [Til They are 
in good agreement with the accurate numerical solutions 
of the previous column. 

To perform an accurate nuclear calculation of the fis- 
sion probability change due to the electron potential it is 
convenient to have an analytical formula for the electron 
energy as a function of the nuclear radius. Therefore, 
we present the result of the fitting of the calculated to- 
tal electron energy of a neutral atom as function of the 
nuclear charge Z and nuclear radius rj^: 



E, 



total 



Eo 



-1) V 10 J 



27 



57(74 

X [0.416 - 0.002(Z - 100)]] , (11) 

where Eq does not depend on nuclear radius: 

Eq = 1.266Z^(7 - l)mc2 [l - 1.52 x 10-'''(Z - 100) 

-2.8 X 10-5(Z- 100)2] . (12) 

These formulas can be used in nuclear calculations of 
the fission process. Note that here we assume Eq. 
for the nuclear charge density p(r) averaged over nuclear 
rotation. This formula contains the fixed parameter D 
and the variable parameter r^r. Nuclear fission calcu- 
lations operate with different variable parameters (e.g. 
the deformation parameters). Relation between r^r and 
these parameters may be easily established by a numeri- 
cal comparison of the mean squared nuclear charge radii. 
It is assumed that nuclear physicists calculate r n for each 
intermediate state of the nucleus during the fission pro- 
cess and then calculate nuclear response to the electron 
potential for a given r^y. To avoid misunderstanding, one 
does not need to calculate the fission lifetime to high ac- 
curacy to make use of the electron effect. One only should 
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calculate the (relative) change of the lifetime when the 
electron potential is added. 

The effect of electrons on the alpha and proton emis- 
sion barriers are much smaller. However, one can use the 
same Eqs.(ll),(12) to estimate them. 

Here we assume that the experiment consists in com- 
parison of the fission probabilities for bare nucleus and 
neutral atom. If necessary, we can provide the result for 
an ion with an arbitrary number of electrons. 

In conclusion we state that we have calculated the 
change of the total electron energy of heavy and super- 



heavy atoms due to a change of nuclear radius and fitted 
the results with simple analytical formula. This formula 
can be used in nuclear calculations to include the effect 
of atomic electrons on the probabilities of nuclear fission. 
Simple estimates based on the parabolic fission barrier 
show that although the effect is small it is probably de- 
tectable in some nuclei. The small value of the effect has 
its advantage: it is enough to consider a linear response 
of a nucleus to the probe - the small rjv-depcndent part 
of the potential produced by the atomic electrons. 
This opens a new way of studying the fission process. 
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